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Abstract 

We study connections between orthogonal polynomials, reproducing kernel functions, and 
polynomials p minimizing Dirichlet-type norms \\pf — l||a for a given function /. For a € [0,1] 
(which includes the Hardy and Dirichlet spaces of the disk) and general /, we show that such 
extremal polynomials are non-vanishing in the closed unit disk. For negative a, the weighted 
Bergman space case, the extremal polynomials are non-vanishing on a disk of strictly smaller 
radius, and zeros can move inside the unit disk. We also explain how dist_D^(l,/ • Vn), where 
Vn is the space of polynomials of degree at most n, can be expressed in terms of quantities 
associated with orthogonal polynomials and kernels, and we discuss methods for computing the 
quantities in question. 


1. Introduction 

The objective of this paper is to study the relationships between certain families of orthogonal 
polynomials and other families of polynomials associated with polynomial subspaces and shift- 
invariant subspaces in Hilbert spaces of functions on the unit disk D = {2 ;GC:|z|< 1}. We 
work in the setting of Dirichlet-type spaces Da, a S K, which consist of all analytic functions 
/ = on the unit disk satisfying 

OO 

( 1 - 1 ) 

Given also g = we have the associated inner product 

OO 

{f, 9)a = + ^Takbk- (1.2) 

fc =0 

We note that Dp C Da when P > a. The spaces Dq, D_i, and Di coincide with the classical 
Hardy space , the Bergman space , and the Dirichlet space D of the disk D respectively. 
These important function spaces are discussed in the textbooks (Hardy), EM (Bergman), 
and [ 10 ] (Dirichlet). One can show that Da are algebras when a > 1, which makes the Dirichlet 
space an intriguing borderline case. Each Da is a reproducing kernel Hilbert space (RKHS): 
for each w G D, there exists an element ka{-, w) G Da, called the reproducing kernel, such that 

f{w) = {f,ka{-,w))a (1.3) 

holds for any / G Da- For instance, when a = — 1, this is the well-known Bergman kernel 
k{z,w) = (1 — wz)~'^. 

Given a function / G Da, we are interested in finding polynomial substitutes for 1//, in the 
following sense. 
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Definition 1. Let / G Da- We say that a polynomial pn of degree at most n G N is an 
optimal approximant of order n to 1// if Pn minimizes \\pf — 1||q among all polynomials p of 
degree at most n. 

It is clear that the polynomials Pn depend on both / and a, but we suppress this dependence 
to lighten notation. Note that for any function / G Da, the optimal approximant Pn (n G N) 
exists and is unique, since p„/ is the orthogonal projection of the function 1 onto the finite 
dimensional subspace f • Vn, where Vn denotes the space of polynomials of degree at most n. 
Note that elements of Da are not always invertible in the space: in general, 1// ^ Da when 
/ G Da- Thus, the problem we are interested in is somewhat different from the usual one of 
polynomial approximation in a Hilbert space of analytic functions. 

Optimal approximants arise in the study of functions / that are cyclic with respect to the 
shift operator S: f i-G- zf. 

Definition 2. A function / G Da is said to be cyclic in Da if the closed subspace generated 
by monomial multiples of /, 

[/] = span{zfc/: fc = 0, 1, 2,...}, 

coincides with Da- 


No function that vanishes in the disk can be cyclic, since elements of [/] inherit the zeros of 
/. The function g = 1 is cyclic in all Da, and if a function / is cyclic in Da, then it is cyclic 
in Dp for all /? < a. If / is a cyclic function, then the optimal approximants to 1// have the 
property 

lbn/-l|U—^ 0: n^oo, 

and the (p„) yield the optimal rate of decay of these norms in terms of the degree n. See 
for more detailed discussions of cyclicity. When a > 1, the algebra setting, cyclicity of / is 
actually equivalent to saying that / is invertible, but there exist smooth functions / that are 
cyclic in Da for a < 1 without having 1// G Da- functions of the form / = (1 — z)^, N G N, 
furnish simple examples. 

In the paper [^, computations with optimal approximants resulted in the determination of 
sharp rates of decay of the norms \\pnf — l||a for certain classes of functions with no zeros in 
the disk but at least one zero on the unit circle T. Thus, the polynomials Pn are useful and we 
deem them worthy of further study. A number of interesting questions arise naturally. For a 
given function /, what are the optimal approximants, and what is the rate of convergence of 
\\pnf — l||a? How are the zeros of the optimal approximants related to these rates, and does 
the location of the zeros of give any clues about whether a function / is cyclic or not? 

In and it was explained that (see Theorem 2.1] for the particular statement 
used here) the coefficients (cfc)^^g of the nth optimal approximant are obtained by solving the 
linear system 

Mc = eo, (1.4) 

with matrix M given by 

(MM)M=0 = (bV,^7)a)M 

and eo = ((1, /),..., (1, z^f)Y' = (/(O), 0). For simple functions /, this system can be solved 
in closed form for all n, leading to explicit expressions for In [^, the authors found the 
optimal approximants to 1/(1 — z) for each a, and plotted their respective zero sets; a plot is 
reproduced in the next section. These plots, as well as the zero sets of optimal approximants 
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for other simple functions / displayed in [^, all had one thing in common: the zeros of the 

polynomials Pn, which we will denote by Z{pn), were all outside the closed unit disk. Might 

this be true for any choice of /, or at least for / non-vanishing in the disk—are optimal 
approximants always zero-free in the disk? 

In this paper we give an answer to this question. For non-negative a, the answer is in the 
affirmative, in a strong sense: optimal approximants are always non-vanishing in the closed 
disk, for essentially any / S Da- 

Theorem A. Let a > 0, let f £ Da have /(O) yf 0, and let (pn) be the optimal 

approximants to 1//. Then Z{pn) flD = 0 for all n. 

For negative a, there is still a closed disk on which no optimal approximant can vanish, but 
this disk is strictly smaller than D. 

Theorem B. Let a <0, let f £ Da have /(O) yf 0, and let (p„) be the optimal 

approximants to 1//. Then Z{pn) H D{0, 2“/^) = 0 for all n. 

We show that the radius 2“/^ cannot be replaced by 1, even if / is assumed to be non¬ 
vanishing, by giving examples of cyclic functions / £ Da, a negative, such that p„(A) = 0 for 
at least one n and at least one A G D \ D{0, 2“/^). 

The proofs of these theorems rely on connections between the Pn, orthogonal polynomials 
in certain weighted spaces determined by the given /, and reproducing kernel functions for 
the polynomial subspaces f -Vn- We also obtain conditions that relate cyclicity of a given 
function / to convergence properties of these orthogonal polynomials and the reproducing 
kernel functions. For example, we show that a function / is cyclic if and only if its associated 
orthogonal polynomials {(pk) have J2k = l/(0)l~^- 

The paper is structured as follows. We begin Section[^by revisiting the optimal approximants 
to l/{\ — z) and by also examining the optimal approximants associated with /jv = (1 — z)^, 

G N; the observations we make in this section motivate much of the further development in 
the paper. 

We point out a connection between the optimal approximants and orthogonal polynomials 
in Section]^ The starting point, given a function / whose optimal approximants we wish to 
study, is to introduce a modified space with inner product {g, h)aj '■= {fg, fh)a', for the Hardy 
and Bergman spaces this amounts to changing Lebesgue measure (on T or D respectively) 
to weighted Lebesgue measure with weight |/p. We study orthogonal polynomial bases for 
the subspace / • Vn and obtain a formula for the optimal approximants in Da in terms of the 
orthogonal polynomials (Proposition |3.1[ ) . For the Hardy space we show that this representation 
implies, via known results concerning zero sets of orthogonal polynomials on the unit circle, 
that the optimal approximants do not have any roots in the closed disk (Theorem |3.2[ ). 

In Section we examine reproducing kernels for the subspaces / • Vn- A relation between 
the reproducing kernel functions and the optimal approximants (see equation ( |4.2[ )) is key, and 
allows us to prove our main result. Theorem |4.2[ 

By combining our results, we can characterize cyclicity of a function in the Da spaces in 
terms of a pointwise (only) convergence property of the sum of absolute values of orthogonal 
polynomials; this is discussed in Section 

Section is devoted to a slightly different idea: the formula ( |1.4[ ) requires the inversion of 
nxn matrices M with (fc, Z)-entry given by = (z^/, z^/). In the case a = 0, the Hardy 
space, multiplication by z* is an isometry, and Mkj = Mk-ip- Hence M is a Toeplitz matrix. 
We use Levinson’s algorithm for inverting Toeplitz matrices to study optimal approximants, 
and we revisit some of the results from the previous sections in the light of this approach. 
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In the last section, Sectionj^ we discuss how zeros of optimal approximants can be computed 
in terms of inner products involving the given function /, and produce examples of functions 
/ S Da, a negative, whose optimal approximants vanish inside the unit disk. 

Some of the results we present here can be readily extended to more general spaces of 
analytic functions, such as Bergman spaces with logarithmically subharmonic weights (see for 
example | 14 [ | 22 [ | 23 [ | 11 | ), but for simplicity, we will concentrate on the Da-spaces as 

defined above. For convenience, we assume /(O) 7 ^ 0 throughout the paper; this simplihes the 
arguments and does not entail any substantial loss of generality. 


2. Motivating examples 

We begin by examining some functions with zeros on T that are cyclic, namely /at = (1 — 
z)^, with N € N\{0}. We present explicit formulas for optimal approximants to l//w and 
investigate their properties, paving the way for further results and conjectures. 


Example 1. For fi = 1 — z, the optimal approximants to 1/fi in Da were found in [^. 
Setting 

Wa(k) = \\z>^\\l = (k + ir, 

the optimal approximants are given by the corresponding Riesz means of nth-order Taylor 
polynomials for 1/(1 — z) = z^. In the series norm for Da that we are considering here, we 
have 


p„(z) = 2^ 1 - ^ ■ 

k=0 \ ^/Wc{j) J 

Using this formula, we can prove the following. 


( 2 . 1 ) 


Proposition 2.1. Let f{z) = 1 — z and let pn denote the optimal approximants tol/f in 
Da- 

(a) The polynomials Pn admit the following representation: 

I ECbAJ*)/ 

(b) The zero set of Pn is given by 

( n+l n+1 'I 

2i{Pn) = <Z Z^/Wa{k) = ^ l/Wa{k) \ . 


( 2 . 2 ) 






(c) In the particular case of the Hardy space (a = 0) the polynomials admit an additional 
representation as follows: 


Pn{z) = 


zn +2 _ (j^j^2)z + n+ l 


(2.3) 


In particular, item (b) tells us that Z{pn) does not intersect D for any n, confirming what 
Figure suggests. Furthermore, an inspection of the formulas reveals that for even n, the 
optimal approximants have no real roots, whereas for odd n, the optimal approximant 
has exactly one real root, which lies on the negative half-axis. 
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Figure 1. Combined zero sets for optimal approximants to 1/(1 — z) in , for n = 0,..., 50. 


Our arguments below are elementary in nature, and clearly limited to this particular / = 
1 — z, and similar functions. Nevertheless, the above observations provided some evidence in 
support of the notion that optimal approximants are zero-free in the unit disk. 


Proof. Parts (a) and (c) can be derived by long division of polynomials: applying Ruffini’s 
rule to the expression (2.1) once yields (2.2 1 and using Ruffini’s rule again on (2.2 1 , gives (2.3). 


Let us verify part (b). Using (|2.2[), we see that Pn{z) can only be zero at singularities of / or at 


points where J2k=o /wk = Since / is entire, and since from (2.1) we know that 


p„(l) yf 0 , we have 


n+1 


n+1 


Z{pn) C := < z ^ 1 : ^ z'^/wk = ^ 1/wk 


k=0 


k=0 


Whenever z € Zn then the numerator in (2.2) is 0 and the denominator is not. Hence Z{pn) = 
Zn- n 


Example 2. We now turn to /at = (1 — z)^, N gN and N >2, which has a multiple root 
at ^ = 1. The optimal approximants to 1//n again admit an explicit representation in the case 
of the Hardy space. If we let B denote the beta function, B{x, y) = — t)y~^dt, then 

the nth-order optimal approximant to 1/(1 — z)^ in is given by 


Pn{z) = ^ 
fc =0 


k + N - l\B{n + N+l,N) 
k ) B{n — k + 1,N) 


zK 


(2.4) 


Once again, in Figure plots suggest that the zeros of the H^-optimal approximants in lie 
outside the closed unit disk for any power N. While this will turn out to be true, we shall see 
in Section that the optimal approximants to 1/(1 — z)^ in Il _2 do vanish in D. 

A proof of Formula 2.4 will be presented in the forthcoming paper |20|, and it seems 


reasonable to suspect that the following holds. 
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Figure 2. Left: Combined zero sets for optimal approximants to 1/(1 — for n = 1,... ,50. 
Right: Combined zero sets for optimal approximants to 1/(1 — z)® for n = 1,..., 50. 


Conjecture 1. The formula (2.4| remains valid for the optimal approximants to l/fa 
when /a = (1 — z:)“, and a G C has positive real part. Note that all such fa are cyclic in the 
Dirichlet space Proposition 13], and hence in all Da for a < 1. 


3. Orthogonal Polynomials 

In order to generalize the observations of the preceding section to arbitrary functions, we 
now turn to a discussion of the relationship between optimal approximants and orthogonal 
polynomials. Fix a G K, and let / G Da, assuming /(O) 7 ^ 0. Consider the space / • Vn, where 
Vn is the space of polynomials of degree at most n. If we let :pkf be an orthonormal basis for 
the space / • Vn, where the degree of ipk is k, then for 0 < fc < n, the functions (pk satisfy 

iVkf-i ^jf}a — ^k,j- 

In other words, we can think of the functions tpk as being orthogonal polynomials in a 
“weighted” Daj space by defining an inner product of two functions tp and in this weighted 
space by 

{‘P,'fp)a,f-= (3.1) 

We let II • \\aj denote the corresponding weighted norm. Without loss of generality we assume 
that each (pk has positive leading coefficient. This choice ensures uniqueness of the functions 
P>k- 


Remark 1. In the case a = 0, where the norm can be expressed in terms of integrals. 


= lim — 
r-^l 2 tt 


\f{re^nVde, 


the space Dqj is simply the weighted Hardy space with dp = \f\'^d9. 

Similarly, when a = —1, the space D-ij is a weighted Bergman space with norm given by 


\\9\\-i,f 


\g{z)\^dp{z), 
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with d/i = l/pdA, where dA denotes normalized area measure. 

For other choices of a, however, equivalent expressions for the norm of 5 S are given in 
terms of the integrals 


l<?(o)P + 


\g'{z)\\l-\z\^f--dA{z), 


and the presence of a derivative means that it is not possible, in general, to write (g, h)aj in 
terms of weighted L^-type inner products in a simple way. 


The optimal approximant minimizes \\pf — l||a over the space of polynomials p G Vn, 
and therefore is the projection of 1 onto / • Vn- Hence, p„/ can be expressed by its Fourier 
coefficients in the basis <pkf as follows: 

n 

{Pnf){z) = '^{l,ipkf)a‘Pk{z)f{z). 
k^O 

Eliminating / from both sides of the expression gives 


Pn{z) = Pkf)aPk{z). 


Notice that by the definition of the inner product (1.2), in all the Da. spaces we have 


{l,g)kf)a = (Pk{0)f{0). 


We have thus proved the following. 


Proposition 3.1. Let a € K and f G Da- For integers fc > 0, let pk be the orthogonal 
polynomials for the weighted space Daj- Let Pn be the optimal approximants to 1//. Then 

n 

Pn{z) = /(O) ^ pk{.b)pk{z)- 
k=0 


Remark 2. Another way to read this expression is as a way to recover the orthogonal 
polynomials pk from the difference between optimal approximants and their values at 0 : 
provided Pn{b) ^ 0 , we have 

, N _ Pn{z) -Pn-liz) 

{^) — ^ . 

<^n(0)/(0) 

We can even recover the modulus of the value at the origin. 


l7>n(0)| 


j PniO) -Pn-l(O) 

W) 


When Of = 0, Proposition |3.1| quickly leads to insights into the nature of zero sets of optimal 
approximants. 


Theorem 3.2. Let / S H^, and let Pn be the optimal approximant to 1//. Then has 
no zeros inside the closed disk. 
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Proof. For / S given, define the positive measure d^{9) = |/(e“)f on the circle, and 

consider the weighted Hardy space of analytic functions g in the disk that satisfy 

r 

\g{e'^)\^dg.{e) < oo. 

Let (/jfc be the orthogonal polynomials for the space iL^(/i), normalized so that the leading 
coefficient Ak of cpk is positive. Now define 


‘Pni.z) ■= Z"(pnil/z), 


(3.2) 


where the polynomial is obtained by taking conjugates of the coefficients of (/?„. Notice that 
if (pn{z) = AnZ"^ + ^ 3 ^^ then (z) = A„ + a^z'^~Z Now it is well-known from the 

theory of orthogonal polynomials (see for example [ 13 [ Chapter 1] or [ 21 [ Chapter 1]) that 


(z) = '^^k{0)‘fk{z). 




A. 


(3.3) 


k=0 


Therefore by Proposition |3.1[ the optimal approximants Pn are multiples of the n—th “reflected” 
orthogonal polynomial: 

Pn{z) = /(0)H„(^*(z). 


Therefore the zeros of are the same as the zeros of 93 *. Moreover it is clear from (3.2) that 


2 ; is a zero of (/?* if and only if 1/z is a zero of pn- Finally, again from the theory of orthogonal 


polynomials, it is well-known that their zeros lie inside the open unit disk (see 13 Chapter 
1 ]), and therefore, the zeros of lie outside the closed unit disk, as desired. □ 

In Sectionwe give a different argument extendable to all values of a. 

In [^, optimal approximants were used to study cyclic vectors, but it is instructive to see 
what happens also in the case when / is not cyclic. 


Example 3 (Blaschke factor in the Hardy space). Let A S D \ {0}, and consider the case 
of a single Blaschke factor 

fx{z) = ^ y ) 

I — Az 

a function that is certainly not cyclic in (or in any Da for that matter). First note that 
\fx{D*)\ = I implies \\pnfx - M\h^ = \\Pn - 1 //a||l 2 , and hence the orthogonal polynomials are 
(fik = z^, k > 0. The optimal approximants are given by 


Pn{z) = '^{l/fx,z'^) 2 Z^. 


k=0 


Note that 1/fx is not analytic in D, but is analytic in C \ D. A calculation shows that 

1//a = A-I-(|Ap — I)z ^-I-• • • 

Therefore, in L^(T), we obtain the coefficients 


{l/fx,z% = 


A k = 0 
0 fc e N+ ■ 


In conclusion, the nth optimal approximant is given by Pn{z) = A for all n, and so is non¬ 
vanishing in the closed disk as guaranteed by Theorem |3.2[ It is not hard to verify 


distnZl, fxVn) = 


\Pn - l/fxl'^dm = 


Y 1 — Az 
A — 


A-. 


dm = I — |A|^. 
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Phrased differently, we have dist// 2 (l, fx ■ Vn) = 1 — /(O)pn(O). In particular, we recover what 
we already know: f\ is only cyclic when A = 1 (and / is interpreted as being constant). 

Both of these observations (non-vanishing of distance formula) will be discussed further 

in the next sections. 


The formula for the n-th reflected orthogonal polynomial expressed in (3.3) relies heavily 
on the fact that / G and that the orthogonal polynomials in this context come from a 
measure defined on the circle. As was explained in Remark no such formula expressing a 
direct relationship between the n-th optimal approximant and the n-th reflected orthogonal 
polynomial holds for measures defined on the disk, and so for Dirichlet spaces Da where a 7 ^ 0, 
such as the Bergman space for example, one must search for different tools. It turns out that 
the language of reproducing kernels is useful in this context. 


4. Reproducing kernels and zeros of optimal approximants 

Let us return to the case of an arbitrary a G M, fix / G Da, let n be a non-negative integer, 
and let (fk be the orthogonal polynomials that form a basis for f ■ Vn, ior 0 < k < n. 

In general, if fc(-,w) is the reproducing kernel function at w in a reproducing kernel Hilbert 
space H, then 

00 

k{z, w) = ^ 'ipk{w)ipk{z) 

for any orthonormal basis {tpk), see [^. Inspecting the relation 
that the function 

n 

Kn{z,w) :='^(pk{w)f{w)ipk{z)f{z) (4.1) 

fc =0 

is the reproducing kernel for the space f ■ Vn- Recall that the reproducing kernel Kn of the 
subspace f ■ Vn C Da is characterized by the property that, for every g G f • Vn, 


3.1 1 now leads to the conclusion 


diw) 


Therefore, by Proposition 
reproducing kernels as follows: 


3.1 


= {g,Kn{-,w))a, WGB. 

the optimal approximants to 1/f are related to these 


Kn{z,0) =pn{z)f{z). 


(4.2) 


One consequence of this fact is the following proposition, whose proof is standard and is 
included for completeness. 


Proposition 4.1. The function Kn{z, 0) / A/Ar„(0,0) is extremal for the problem of finding 

sup{|ff(0)| :gGf-Vn, ||g|U < 1}: 
and thus the supremum is equal to x/ Kn{0,0 ). 

Proof First note that ||iL„(-, 0)||^ = (Ar„(-, 0), Ar„(-, 0))^ = Ar„(0,0), by the reproducing 
property of Kn{-, 0). Now let g be any function in f ■ Vn such that ||g||a < 1. Then 

|<?(0)| = \{g,Kn{;0))a\ < \\g\U ll^n(-,0)|U < \/^n(0,0). 

Choosing g{z) = Kn{z, 0)/-\/Ar„(0,0) gives that ||g||a = 1 and g{0) = x/Kni0,0), and thus g is 
a solution to the extremal problem stated in the proposition, as required. □ 
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Expressing the optimal approximants in terms of these kernels allows us to prove our main 
result concerning zeros of optimal approximants. 

Theorem 4.2. Let a G M., let f G Da have /(O) 7 ^ 0, and letpn be the optimal approximant 
to 1/f of degree n. Then 

- if a > 0 , all the zeros of the optimal approximants lie outside the closed unit disk; 

- if a < 0, the zeros lie outside the closed disk D{0, 2“/^). 


It is clear that the kernels Kn vanish at all the zeros of /. Borrowing terminology from 
Bergman space theory, we say that any A G C such that iG„(A,0) = 0 but /(A) 7^ 0 is an 
extraneous zero. Theorem 4.2 can then be rephrased by saying that the reproducing kernels 
Kn{-,0) have no extraneous zeros in D when a > 0, and no extraneous zeros in 11(0, 2“/^) when 
a < 0. 

Proof Let kn{z) := Kn{z,0) =p„(z)/(z) be the reproducing kernel at 0 for the space / • 
Vn, and suppose A is an extraneous zero of kn- Then 

kn{z) = {z- X)q{z)f{z), 

where g is a polynomial of degree at most n — 1. Therefore 

= zq{z)f{z) - kn{z). 

Notice that since zq G f ■ Vn and vanishes at 0, while kn reproduces at 0, 

0 = {zqf){0) = {zqf,kn)a, 
the two functions zqf and kn are orthogonal. It follows that 

mqf\\l = \\zqf\\l + \\kn\\l (4.3) 

For any function F(z) = ^nZ^ in Da, we have 

00 

\\F\\l = Y.(^ + ir\an\^ 

while 

00 

\\zF\\l = Y,{n + 2r 

It is clear that 


n—0 


n—0 


E 

n—0 


n + 2 
n + 1 


{n + iy 


Tl -\~ 2 

1 < - < 2 . 

n + 1 

Hence, if a > 0, then || 2 :T||^ > IjTlI^, while if a < 0, we obtain || 2 T||q > 2“||T||^. Applying 
these estimates to F = qf in ( |4.3[ ), we obtain, for a > 0, that 

{\M^-i)hf\\l>\\kn\\l>0, 

which implies that |A| > 1, as claimed. For a < 0, it follows that 

{\X\^-2-)\\qf\\l>\\kn\\l>0, 

which implies that |A| > 2“/^, as desired. □ 


A few remarks are in order. 






ORTHOGONAL POLYNOMIALS, KERNELS, AND ZEROS 


Page 11 of[^ 


Remark 3. Note that in the case of the Bergman space another way to see the relationship 
between the norm of a function and the norm of its multiplication by 2 is to recall that 
the function \ fl z is the so-called “contractive divisor” at 0 , and thus is an expansive 
multiplier (see 14 or [^). Therefore one has Hv^zFU-i > ||F||_i, which is equivalent to 
the desired inequality. The same remark applies to in the range a G [—2,0]. Moreover, it is 
straightforward to show (see, e.g., 0 ) that when n —>■ oo, kn{z) —>■ k{z, 0 ), where k{z, w) is the 
reproducing kernel in the weighted space Daj, when / is sufficiently nice up to the boundary 
of the disk. As is known (i). fc(z, 0) has no extraneous zeros. Thus, for a G [—2,0), the zeros 
of kn are all eventually “pushed out” of the unit disk when n —> oo. 


Remark 4. The proof of Theorem 4.2 is similar to a well-known proof (due to Landau, 
according to 21 ) about the location of the the zeros of orthogonal polynomials in a fairly 
general setting. For example, suppose /it is any measure on the unit disk D and let be the 
orthogonal polynomial of degree n with respect to /r, normalized for instance by requiring its 
leading coefficient to be positive. Then 


(Pn{z)(Pm{z)d^,{z) = 0 
B 


\i n ^ m, and so ipn is orthogonal to any polynomial of degree strictly less than n. Now if 
^n{^) = 0, we can write (pn{z) = (A — z)q(z), where g is a polynomial of degree n — 1. Then 
zg(z) = \q{z) — q)ri{z), and therefore 

ll^gf = |An|gf+ ||^„f. 


Since 2 ; G D, ||zg|| < ||g||, and therefore we obtain that 

(l-IAr)llgf >0, 

which implies that |A| < 1. In fact, one could refine the estimate further based on the support 
of fj,, for instance if /r were an atomic measure, since 

jjzgf < max{|z| : z G supp(/i)} • ||gf, 


one would obtain that |A| < max{|z| : z G supp(/x)}. 


Remark 5. We do not know whether the radius 2“^^ is optimal, that is, whether there 
are examples of optimal approximants, associated with functions in with a negative, that 
vanish at points A with modulus arbitrarily close to 2“/^. 

However, in Section]^ we present examples of functions that lead to extraneous zeros located 
at A « 0.88 when a = —2, and at A » 0.9 8 when a = —1. Hence D{0, 2“/^) cannot be replaced 
by D in the second part of Theorem 


4.2 


5. Conditions for cyclicity 

We can use the reproducing kernels to give equivalent criteria for the cyclicity of / in terms 
of the pointwise convergence of kernels at a single point, the origin. 

Theorem 5.1. Let f G Da, /(O) 7 ^ 0, and let pn be the optimal approximant to 1/f of 
degree n. Let (pk be the orthogonal polynomials for the weighted space Daj- The following are 
equivalent. 

(i) / is cyclic. 

(ii) p„( 0 ) converges to l//( 0 ) as n —)• 00 . 
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(iii) j:Zo\Mo)\^ = y\fm^- 


In the next Section we will extend this theorem in the case when a = 0, by including an 
additional equivalent condition. 


Proof. We first record some observations. By (4.2) we find: 


/ is cyclic •, 0) - 1| 


0 . 


(5.1) 

(5.2) 

(5.3) 

Now we show (1) ^ (2) ^ (3) ^ (1). 

Assume (1). Then (2) follows from pointwise convergence of to 1//. By (4.2) we obtain 
Kn{0,0) —t 1, which implies item (3) by virtue of (5.2). 


From Equation (4.1), we obtain 

n 

\\Kn{ -, 0)112 = K^{0,0) = |/( 0)|2 ^ \M0)\". 

k^O 

Now, for any function h in we have the orthogonal decomposition of the norm as 

\\h\\i = \hm^ + \\h-hmi 


It remains to argue that (1) follows from (3). Assume that (3) holds. Then by (5.2) as n —>■ oo: 


Kn{0, 0) —1, and 

l^n(-,0)||^l. 


(5.4) 

(5.5) 


Equation (5.3) applied to h = Ar„( • ,0) yields 


\Kni-,0)\\l = |if„(0,0)|2 + \\K^{-,0) - if„(0,0)||2. 


In view of (5.4) and (5.5) we learn 

||iF„(-,0)-i^„(0,0)|U^0. 

Since monomials are orthogonal in 

\\Kn{ • , 0) - 1||2 = ||iF„( . , 0) - X„(0, 0)||2 + |iA„(0, 0) - 1|2 


(5.6) 


and this, again together with (5.4), informs us that ||Ar„( • ,0) — 
from (5.1). 


0. Item (1) now follows 
□ 


A further equivalent criterion can be formulated by relating the distance distD^(l,/ • Vn) 
to the values of pn at the origin, as in Example This is actually part of a general statement 
contained in a classical result of Gram, which we phrase here in our terminology, although 
it applies in any Hilbert space, and for distances involving more general finite-dimensional 
subspaces. 

For / € Da fixed, consider the matrix 

Mn = = {{z^ f,Z f)a)jj,^Q 

and denote its lower right n-dimensional minor by Mn- 

Lemma 5.2 (Gram’s Lemma). Let f G Da- Then = dist£)^(l,/ • Vn) satisfies 

dl = l-pn{0)f{0) (5.7) 

where Pn is the nth optimal approximant to 1//. Moreover, p„(0) is given as 

Pn{0) = /(O) det Mn /det Mn- 


(5.8) 
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Proof. First, notice that, since the orthogonal projection of 1 onto Vnf of 1 is Pnf, we have 

dl= {Pnf -f,Pnf -f)a- ( 5 - 9 ) 


Since Pnf — f is orthogonal to all functions of the form qf where q is a polynomial of degree 
less or equal to n, we obtain that 


dl = {Pnf - 1, -l)a = 1 - Pn(0)/(0). 


(5.10) 


The matrix satisfies M„c = b where c are the coefficients of and bk = (1,z^/)q. These 
equations together with (5.10) form a system of equations with unknowns and c. Using 
Cramer’s rule to solve for gives the determinant identity (5.8). □ 


Gathering everything we have obtained so far, and including Gram’s Lemma, we obtain the 
following. 


Corollary 5.3. Let / g Da satisfy /(O) ^ 0 and let pn, and Kn be as above. Then 
the following quantities are all equal: 

(a) dist^^(l,/• T’n) 

(b) \\Pnf-nl 

(c) l-p„(0)/(0) 

(d) l-((M„)-i)o,o|/(0)P 

(e) l-ELo\MO)?\f{0)? 

(f) l-d^„(0,0) 

If, moreover a = 0 and the degree of Pn is equal to n, then all the above are also equal to 


Hence, / is cyclic if any (hence all) of these quantities tend to zero with n. Since the distance 
in (a) above is always a number between 0 and 1 and converges (as n goes to oo) to the distance 
from 1 to [/], and the numbers in item (e) are non-increasing, all the other quantities converge 
in the interval [0,1]. In particular, a function / is cyclic if and only if the kernel of the invariant 
subspace generated by /, K[f], satifies df[j](0,0) = 1. 


Remark 6 (A formula of McCarthy). 

We point out a connection with an observation of McCarthy, see [ 18 |[ Theorem 3.4]. Under 
the assumption that / G H°° is cyclic in the Bergman space (a = — 1), he provides a closed 
formula for the reproducing kernel K of the closure of the polynomials with respect to jj • ||-i,/. 
His result generalizes to Da with a < 0, and yields 


K{z, w) 


1 1 
7H/(z) (i-wz)i-“’ 


z, ru G D. 


This is in effect a rescaling of the reproducing kernel of Da] see Ghapter 2.6] for a discussion 
of this notion. 


6 . Toeplitz matrices and Levinson algorithm 


In view of (1.4) and Gorollary 5.3 it is of interest to consider different algorithms for inverting 
the matrices M„. 

Multiplication by is an isometry on Therefore, in the case of the Hardy space, the 
matrices M = (M^j) appearing in the determination of the optimal approximant have the 
property that M}^ i = In other words, the entries M]^ i only depend on the distance 
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to the diagonal. A matrix with this property is called a Toeplitz matrix. We can use this 
structure of the matrices to extend our results in Theorem |5.1[ A number of algorithms have 
been developed specifically for inverting Toeplitz matrices. In 16 , p. 7-13] several methods 
are mentioned, based either on Levinson’s 17 or Schur’s 19 algorithms. 


Theorem 6.1. Let f G be such that /(O) ^ 0, and denote Ck,r. 
the n-th optimal approximant to 1//. Then f is outer if and only if 


the k-th coefficient of 


n 

n—O 


1 - 


|Cn+l, 

|co, 


n+l| 


n+l| 


= /(0)/||/ll 


Remark 7. Notice that, in the notation for orthogonal polynomials used in the previous 
sections, c„^„/co,n = <p„(0)/(^„(n), and this quotient is also the product of the numbers 
where Zk varies over all zeros of p„, or alternatively, the product of the zeros of the 
orthogonal polynomial ipn- In particular, in the Hardy space, cyclicity can be characterized 
exclusively in terms of the zeros of optimal approximants or in terms of those of orthogonal 
polynomials. Theorem |6.I| is, in a sense, a qualitative optimal approximant version of the known 
characterization of outer functions as those satisfying log |/(0)| = ^ log \f\d9. It would be of 
great interest to know whether a version of Theorem |6.1| also holds in other spaces. 


Proof 

Without loss of generality we assume /(O) = 1 (otherwise divide / by /(O)). As was explained 
in the introduction, the coefficients c = (cq, Ci,..., c„) of the optimal approximant of order n 
are given by the linear system 

Mc=eo where eo = (l,d) and Mk^i = {z^f,z^f). 

By virtue of the existence and uniqueness of the minimization problem, the matrix M is 
invertible. Our objective is to obtain the coefficients by taking 

c = M~^eo. 


Now we will use the fact that M is a Toeplitz matrix, and apply the Levinson algorithm. As 
our matrix is in fact Hermitian, we can apply a slightly simplified version of this procedure. The 
algorithm is based on the fact that all information of the matrix is contained in two columns 
(when M is Hermitian, in one column). 

The solution is as follows: If {ck.n}k=o coefficients of the nth-degree optimal 

approximant, then the coefficients {ck,n+i}^t,li of the optimal approximant of degree n-I-1 
can be obtained from those previous coefficients: 


where 


Ofc,n-t-l — iT^ 19 (Ofc.n L.fiCji.\-i — k , 

1 - Tnr 

n 

r„ = ^c„_fc,„(z'=+V,/)- 

k—0 


( 6 . 1 ) 


Since /(O) 
obtain 


1, the numbers Co,n are always real. From the expression above, we can then 


-r„ 


1- |r. 


, eo,n — TnCQ^ri+l- 


Cn+l,n+l — 


( 6 . 2 ) 
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Finally, this gives us 


r„ = - 


Cn+l,n+l 

C0,n+1 


From (6.1) we can recursively recover the value ofp„(0): 

Co,0 

Pn+l(0) - |p^|2 


(6.3) 


(6.4) 


nLo(i-ir.P)’ 

where r„ is defined as in ( |6.3[ ). The value of co,o can be recovered from the corresponding 
equation ||/||^co,o = 1 and thus, (6.4) becomes 

Pn+l(0) = YrLii^’fT" IT^ |2\ ■ (®-^) 

li/r nfe=o(l “ Tfel ) 

Being outer is equivalent to cyclicity in Hardy space by the classical theorem of Beurling. 
By Theorem 5.1 this will happen if and only if p„(0) tends to 1 as n tends to infinity, which 
happens if and only if 


n(i-|rfe|2) = 1/11/11^ 


fc =0 


Using (6.3) again to translate the value of Ffc, we obtain the desired result. 


□ 


Example 4 (Optimal approximants to 1/(1 — z) revisited). We illustrate how the Levinson 
algorithm can be exploited for our purposes by using it to re-derive the optimal approximants 
for the basic example / = 1 — z. In this case, the main advantage is that F^ is very simple: 

Pn — ^ G M. 

So to verify that Cesaro polynomials, which correspond to the choice Ck^n = 1 — (fc + l)/(n + 
2), are optimal, we just need to check that they satisfy the recursive formula and the initial 
condition for degree 0. 

Hence, we want to show that 

Cfc,n+1 — Z I TZ (Cfc n 4“ ^n^n^n+l — k,n) • 

1 - |Cn,nr 

Evaluating both sides of the formula reduces our task to checking that 

n + 2 — k (n + 2)(n + 1 — k) -\-k 
n + 3 (n + 3)(n + 1) 

Multiplying both sides by (n + 3)(ri + 1), we obtain that 

(n + l)(n + 2 — fc) = + 3n + 2 — nA: — fc = (n + 2)(n + 1 — k) + k, 

and therefore the proposed polynomials are optimal as claimed. 


Remark 8. A particularly simple case is that of / an inner function in the Hardy space. 
Then, in the notation of the proof of the previous Theorem, Ffc = 0 for all k > 0, and the 
polynomials do not change with k. This means that either the optimal norm converges to 0 with 
Pn all being equal to a constant, or it does not converge. In fact, when f = IF, where I is inner 
and F, outer, the elements of the system (1.4) depend exclusively on F. That is, for / S H'^, 
the optimal approximants depend only on the outer part of /. Another immediate consequence 
is that for any outer nonconstant function /, there is some t > 1 such that (z*/, /) 0. In 

other words, inner functions are characterized by having optimal approximants of all degrees 
equal to a constant. 
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7. Extraneous zeros 

We now return to the zero sets of optimal approximants, with a view towards determining 
the location of zeros analytically. Our first result states that the roots of Pn can be expressed 
in terms of certain inner products. 


Lemma 7.1. Let / S Da have /(O) ^ 0,and let Z{pn) = {zi ,..., z„}. 

Then are given by the unique-up-to-permutation solution to the system of 

equations 


— 


- Zj) ^ 

{f^]=i,jym(z - Zj),zfU^^^^^^^{z - Zj)) 


(7.1) 


for m = 1 ,..., n. 

In particular, the zero ofpi, the first order approximant, is given by 

WzfWl 


zi = 


if, zf). 


(7.2) 


If some zero Zm is repeated, the solution is still unique but we count multiplicity. Note also 
that if / is a cyclic function but not a constant or a rational function, then there is an infinite 
subsequence {nk}k&-i such that deg(pnfc) = rik; if not, \\pnf — l||a cannot tend to 0 as n —>■ oo. 


Proof. The first-order approximant pi = cq + Ciz to 1// is obtained by solving the system 
of equations 

ll/ll"co + (^/,/)„ci=7M (7.3) 

and 

{f,zf}aCo + \\zf\\lci =0. (7.4) 

Suppose Cl ^ 1 (otherwise interpret zq as being equal to infinity). Then p{zi) = 0 is equivalent 
to 


and by (7.4) then, we obtain 


zi = -cqIci- 

WzfWl 

(/. Zf)a 


Zl = 


Next, we note that (7.2) can be expressed as the orthogonality condition 

{{z- zi)f,zf) = 0. 


(7.5) 


To prove the lemma for the optimal approximant of 1// of any degree, it is enough to apply 
Equation (7.5) to a function g that is the product of / with a polynomial of degree n — 1: 

{fP%^{z - z,),zfll]^,^^^,iz - Zj)) = 0 . 


If the optimal polynomial to invert / has n zeros, each of them has to satisfy a corresponding 
orthogonality condition for a different g. Moreover, the fact that we multiply the polynomial by 
a constant does not affect the orthogonality condition, and the zeros are determined exactly by 
those orthogonality conditions. Since the polynomials are unique, the zeros are also uniquely 
determined. □ 


We shall now use Lemma 7.1 to show that optimal approximants to 1// have zeros in D for 
judiciously chosen /, or in other words, that the associated kernels have extraneous zeros. 
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Figure 3. Zeros of the optimal approximants pi (black), p 2 (green), and ps (gold) associated 
with (1 + z)^ in the weighted Bergman space I?_ 2 - 


We present two families of examples, one that is completely elementary, and one that requires 
more work but has the advantage of producing extraneous zeros in the disk for the classical 
Bergman space. 


Example 5 (Extraneous zeros in weighted Bergman spaces). We begin by treating the 
spaces Da with a < — 1. An equivalent norm for Da is given by the integral 




B 

and so the spaces coincide (as sets) with the standard weighted Bergman spaces discussed 
in (Hl§ and also studied in 


22 


_ _among other references. 

We return to the functions /at = (1 + z)^ and set N = 3. A direct computation shows that 
the first optimal approximant to 1/(1 + in D -2 vanishes at 

11^(1+ 11-2 741 


° {{1 + zY, z{\ + zY) -2 755 

a point inside the disk. By differentiating the function 

11 ^( 1 + 11 ^ 


= 0.981. 


((1 + Z)3,Z(1 + Z)3), 


with respect to a, we see that zq is increasing on the interval (—oo, —2], and hence the optimal 
approximant pi to 1/(1 + z)^ has a zero in D, for any Da with a < —2. In fact, by choosing 
N = N{a) large enough we can produce an extraneous zero also for the range —2 < a < 1. We 
omit the details. 

Straight-forward linear algebra computations produce the first few optimal approximants to 

1/(1+ zY: 


741 / 775 \ 

1694 V ~ ^ ’ 


P2 


961 / 1571_ 1032 

1(338 ~ ) ' 


Pi = 
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Figure 4. Zeros of the optimal approximants pi (black), p 2 (green), and ps (gold) associated 
(1 + z)^ in the unweighted Bergman space H-i. 


and 


571 / 3427 1182 , 



It can be checked that the zero sets Z(p„), n = 1,2,3, are all contained in the unit disk; see 
Figure]^ 

The second source of examples is the family of functions 


/.= 


1 + z 
( 1 - 2 )"’ 


r] > 0. 


We have 


/i 


1 + z 
1 - z 


1 + 25] z '= 

fc=i 


and we see that /i G 71-2■ Moreover, /i is cyclic as a product of the cyclic multiplier 1 + z and 
the function 1/(1 — z), which is cyclic in for allp < 1, and hence also in I?_ 2 . Using Euler’s 
formula = 7r^/6 to compute ||z/i ||?.2 and (/i,z/i)_ 2 , we find that the first-order 

optimal approximant to l//i in D _2 vanishes at 


87r2 - 57 
87r2 - 54 


0.879... 


Example 6 (Extraneous zeros in the Bergman space). It can be checked, again by 
hand, that the zeros of the first few optimal approximants to 1/(1+ z)^ in the Bergman 
space D_i are in the complement of the unit disk, see Figure]^ In fact, one can show that 
I|2/a||^i/(/a,^/a)-i > 1 for all NgN. 

However, this is not always the case! Before presenting a specific example, let us give a 
heuristic explanation for why the zeros of optimal approximants may move inside D for Bergman 
type spaces. Let / be a cyclic function in the Bergman space 7?_i, say, and define g(z) = 
z/(z), assuming the normalization || 5 |j_i = 1. Then by (7.2), we need to find / such that 











ORTHOGONAL POLYNOMIALS, KERNELS, AND ZEROS 


Page 19 of[^ 


\zi\= < 1 : or equivalently, 


-\g\'^dA 

z 


> 1 . 


Letting h be defined by 5 = y^zh, this equation becomes 

z\h\'^dA > 1/2, 


(7.6) 


where, since z is a contractive divisor and hence an expansive multiplier in all Bergman 
spaces with logarithmically subharmonic weight (cf. mi) , we have that ||/i|l-i < 1- In other 
words, we are looking for h such that the measure dg := j/ipdA has total mass at most 1 but 
has center of mass close enough to 1 to ensure that (7.6) holds. Thus if we are able to choose 
h so that fi is concentrated in the circular segment 5" := {z S D : Re(z) > 1 — e} for small e, 
and say h is symmetric with respect to the x-axis, then the center of mass of /i will be real and 


close to 1, so inequality (7.6) will be satisfied. Starting with f{z) = for 0 < /? < 1 but 

sufficiently close to 1 , all the requirements will be fulfilled and the zero of the first approximant 
will move inside D. 

The following example is adjusted from the above idea to make the calculations come out in 
essentially closed form. Specifically, let us consider the function = (1 + z)/{l — z)'^/^. We 
note that is cyclic in U_i since 1 + z is a cyclic multiplier, and V (1 ~ is cyclic in 
the Hardy space which is contained in the Bergman space (see [9]). 

By the binomial theorem, we have frj = l + ^k{il)z^, with 


Ckiv) = (- 1 )^ 


k 


-v 

k-l 


Using a computer algebra system, such as Mathematical one checks that 

2 


Aiv) = Y. 


-V 


fc=i 


k + 2 


„ .,T(2-2r]) 1 

“ ^ ^[r(3-7?)]2 " 2’ 


B{r]) = Y ^ , o ^ = I 3 F 2 { 3 , r], 77 ; 1, 4; 1), 


and 


civ) = -Y 


k=l 


-V 

k-l 


k + 2 


fc=i 


k + 2 


1 [ r(2-2r;) 

1-77 \T{1 - fj)Ti2r,) 


- 3-F2(2,77 - 1,77;1,3;1) I . 


Here, 3 F 2 denotes the generalized hypergeometric function. Evaluating at 77 = 4/5, expressing 
everything in terms of gamma functions and repeatedly using the functional equation r(a; + 
1 ) = a;r(a:), we find that 

■1(4/5) = g 


and 


Upon combining, we obtain 


25 [r(ll/5)]2’ 265 [r(16/5)]2’ 

2464 r(2/5) 


C(4/5) = 


625 [r(16/5)]2‘ 


N/usIlii = 


2 _ 2142 r(2/5) 


125 [r(16/5)]2' 
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A similar analysis applies to (/ 4 / 5 , 2 :/ 4 / 5 )-i, and we find that 

r(7/5) 


(/ 4/51 zf4i5)-i — 9 

After simplifying the resulting ratio, we obtain 


[r(ii/5)]2' 


Zi = 


2 / 4 / 5 II -1 


119 


(/ 4 / 512 / 4 / 5 )-! 121 

and so pi has a zero in the unit disk, as claimed. 


= 0.983..., 


It is possible that, with additional work, one could use to exhibit extraneous zeros also 
for Da in the range — 1 < a < 0, but this seems more technically challenging. 


Remark 9. The failure of Bergman space analogs of results for Hardy and Dirichlet spaces 
is a common occurrence. One example of this phenomenon that seems relevant is the existence 
of non-cyclic invertible functions in the Bergman space that was discovered in [3. This is in 
contrast to and the Dirichlet space, where invertibility implies cyclicity. In^, as in our 
Example]^ the source of unexpected bad behavior is not, as one might predict, a “large” set on 
the boundary where the function vanishes, but rather the presence of regions of rapid growth 
of the function. 

Another example, close in spirit to the previous example, of how Hardy and Bergman 
spaces are different can be found in [ 12 | . There, it is shown that while eigenfunctions of a 
certain restriction operator acting on ID never vanish on the unit circle, eigenfunctions of the 
corresponding operator on the Bergman space may indeed vanish on T. We thank Harold S. 
Shapiro for pointing out this reference to us. 

Viewed in a different light, it is perhaps somewhat surprising that there are extraneous 
zeros inside the disk in the case of the unweighted Bergman space. An important step in 
the construction of contractive divisors for Bergman spaces (see 


14 


) is to rule out 

extraneous zeros of a certain extremal function. This can be done for the Bergman space, and 
more generally for the weighted spaces Da for —2 < a < 0, but extraneous zeros do appear 
when a < — 2, see [ 15 |. In our case, zeros in the disk are present already for a = — 1. 


Remark 10. Since all of our examples are cyclic vectors, the associated reproducing kernels 
have to converge, as n —>■ 00 , to the reproducing kernels of the respective Da- These latter 
kernels are zero-free, and hence the zeros of have to leave every closed subset of the unit 
disk eventually. See for details. It does not seem easy to determine how fast this happens, 
or whether there is any monotonicity involved: in principle it could happen that some Pn is 
zero-free, while some subsequent Pn' again vanishes inside D. It is known, see Section 3], that 
monotonicity does not hold for zeros of Taylor polynomials associated with outer functions. 
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